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Abstract 



Toral (2002) considered an ensemble of N > 2 players. In game B a 
player is randomly selected to play Parrondo's original capital-dependent 
game. In game A' two players are randomly selected without replacement, 
and the first transfers one unit of capital to the second. Game A' is fair 
(with respect to total capital), game B is losing (or fair), and the random 
mixture "fA' + (1 — 'y)B is winning, as was demonstrated by Toral for 
7=1/2 using computer simulation. We prove this, establishing a strong 
law of large numbers and a central limit theorem for the sequence of 
profits of the ensemble of players for each 7 £ (0, 1). We do the same 
for the nonrandom pattern of games {A'Y B" for all integers r,s > 1. 
An unexpected relationship between the random-mixture case and the 
nonrandom-pattern case occurs in the limit as A'^ ^> 00. 

AMS 2000 subject classification: Primary 60J10; secondary 60F05. 

Key words and phrases: Parrondo's capital-dependent games, Markov 
chain, stationary distribution, fundamental matrix, strong law of large 
numbers, central limit theorem. 



*This work was partially supported by a grant from the Simons Foundation (209632 to 
S.N.E.). It was also supported by the Korean Federation of Science and Technology Societies 
Grant funded by the Korean Government (MEST, Basic Research Promotion Fund). 

tThis work was supported by the Basic Science Research Program through the National 
Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and 
Technology (2011-0005982). 



1 Introduction 



The original Parrondo (1996) games can be described in terms of probabilities 
p ■.— 1/2 — £ and 

Po:=Yg-£, pi=p2-= ^-e, (1) 

where e > is a small bias parameter (less than 1/10, of course). In game A, 
the player tosses a p-coin (i.e., p is the probability of heads). In game B, if the 
player's current capital is congruent to j (mod 3), he tosses a pj-coin. (Assume 
initial capital for simplicity.) In both games, the player wins one unit with 
heads and loses one unit with tails. 

It can be shown that games A and B are both losing games (asymptotically), 
regardless of e, whereas the random mixture {1/2){A + B) (toss a fair coin 
to determine which game to play) is a winning game for e sufficiently small. 
Furthermore, certain nonrandom patterns, including AAB, ABB, and AABB 
but excluding AB, are winning as well, again for e sufficiently small. These are 
the original examples of Parrondo 's paradox. 

It has been suggested that game A acts as "noise" to break up the losing 
cycles of game B played alone (Kay and Johnson 2003). Toral (2002) proposed 
a stochastic model in which a different type of noise appears to have a similar 
effect. The model assumes an ensemble of iV > 2 players and replaces the noise 
effect of Parrondo's game A by a redistribution of capital among the players. A 
player i is selected at random to play. With probability 1 /2 he can either play 
Parrondo's game B or game A' consisting in that player giving away one unit of 
his capital to a randomly selected (without replacement) player j. Notice that 
this new game A' is fair since it does not modify the total amount of capital, it 
simply redistributes it randomly among the players. 

Toral showed by computer simulation that the Parrondo effect is present in 
his games. Our aim here is to prove this, establishing a strong law of large 
numbers and a central limit theorem for the sequence of profits of the ensem- 
ble of N players. For this we apply results of Ethier and Lee (2009), but the 
application is not straightforward. For example, the formulas for the mean and 
variance parameters in the central limit theorem depend on the unique station- 
ary distribution of the underlying Markov chain as well as on its fundamental 
matrix, both of which are too complicated to derive explicitly except for small 
N. Nevertheless, we can evaluate the mean and variance parameters for all N. 

We generalize ([T]) to the parameterization of Ethier and Lee (2009): 

1 + p"^ 1 + p 

where p > (eq. (P) is the special case p = 1/3). The bias parameter is not 
important, so we take £ = in most of what follows, which makes game B fair 
(asymptotically) . 

Let us summarize our results. Just as it is conventional in the literature to 
denote the nonrandom pattern {A'YB'' by [r,s], we will introduce the (slightly 
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redundant) notation (7, 1 — 7) for the random mixture 7A' + (1 — 'y)B. We 
establish a strong law of large numbers (SLLN) and a central limit theorem 
(CLT) for the sequence of profits of the ensemble of N players in both set- 
tings (random mixture and nonrandom pattern) . Wc provide a formula for the 
random-mixture mean ^J■\^^\_^y, which does not depend on A^, as a function of 

7 € (0, 1) and p > 0. The nonrandom-pattern mean /^[^jj does depend on N and 
is rather more complicated; we provide a formula, as a function of iV > 2 and 
p > 0, only for small r,s > 1 but we determine its sign for all r, s > 1, > 2, 
and p > 0, thereby establishing necessary and sufficient conditions for the Par- 
rondo effect to be present. Finally we show that the random-mixture case and 
the nonrandom-pattern case are connected by the unexpected relationship 

p'}^}, ^ , lim uf^^l, r,s> 1, N > 2, p> 0, (3) 

and a simple formula for this common value is provided. To put this in perspec- 
tive, the corresponding identity for one-player Parrondo games fails in all but 
one case (r = 2, s = 1). 

The variance parameter is considerably more complicated, so we assume 
that p = 1/3 (i.e., ([l} holds with e = 0) and 7 = 1/2, obtaining a formula 
for (cl^^j 1/2))^ ^' function of > 2. We do the same for (cr'^^^^)'^ for p = 
1/3 and small r,s > 1. It turns out that the analogue of fails for the 
variances. However, a different notion of variance, the expected sample variance 
of the individual players' capitals, which was considered by Toral (2002), does 
apparently satisfy a relationship nearly analogous to ([3|). We can confirm this 
only in special cases, so it remains a conjecture. 

Toral (2002) also studied a model in which the capital-dependent game is 
replaced by a history-dependent game. It seems likely that most of the results 
of this paper can be extended to that setting, with the probable exception 
of Theorem [7] below. Finally, Toral considered a model with redistribution 
of wealth from richer to poorer neighbors. That model is considerably more 
difficult to analyze than either of the others, and no rigorous results for it are 
known. 



2 Mean profit for random mixtures 

There are two natural ways to define the model. The simplest is to describe the 
state of the system by an A^-dimensional vector x = (xi,X2, . . . ,xn) in which 
Xi denotes the capital (mod 3) of player i. An alternative approach (adopted by 
Ethier 2007), which makes the state space smaller but the one-step transition 
probabilities more complicated, is to describe the state of the system, when it is 
in state x according to the previous description, by (rip, rii, 712), where uq (resp., 
ni, 71,2) is the number of Os (resp.. Is, 2s) among xi, 2:2, . . . , xjy. Using the first 
approach, the state space is 

Eat := {x = {xi,X2,..., xn) : G {0, 1, 2} for i = 1, . . . , iV} = {0, 1, 2}^, 
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while using the second approach, the state space is 

Sat := {(no, ni, 712) € : no + ni + n2 = N}. 

We note that |EAr| = 3^ and ^ {^+^). 

The one-step transition probabihtics using the first approach depend on three 
probabilities po,pi-,P2- If only game B is played, then the one-step transition 
probabilities have the simple form 

P^^\x y) J -^"^^^^ ii = Xi + 1 (mod 3) and yj = for ah j ^ i 
^ ' I N~^qx- if j/i — Xi — 1 (mod 3) and yj — Xj for all j ^ i 

for i ~ l,2,...,iV, where qx '■— ^ — Px for x = 0,1,2, and Pl^\x,y) ~ 
otherwise. We adopt the parameterization ^ with e = 0. 

If only game A' is played, then the one-step transition matrix is symmetric 
and of the form 

P[':'\x,y) := [N{N-l)]-' 

if, for some i,j S {1,2, . . . , N} with i ^ j, we have yi = Xi — 1 (mod 3), 
j/j = Xj + 1 (mod 3), and yk = x^ for all k ^ Finally, if the two games 
are mixed, that is, game A' is played with probability 7 € (0, 1) and game B is 
played with probability 1 — 7, then our one-step transition matrix has the form 

The one-step transition probabilities using the second approach also depend 
on the three probabilities pQ,pi,p2 and are best summarized in the form of a 
table. See Table [1] which is essentially from Ethier (2007). 

That the two approaches to the model are equivalent, at least in the sta- 
tionary setting, is a consequence of the following simple lemma, which is easily 
seen to be applicable to P^^^ and P^^^_^y 

We first need some notation. Given a finite set E and an integer N > 2, 
put := E X ■ ■ ■ X E. Given a permutation a of {1, 2, . . . , A^} and x = 
{xi, . . .,xn) e -B^, write x„ := (x^(i), . . .,x„^n))- 

Lemma 1. Let E be a finite set, fix N > 2, let P be the one-step transition 
matrix for an irreducible Markov chain in the product space E^ , and let tt be 
its unique stationary distribution. If, for every permutation a 0/ {1, 2, ... , N}, 

P{xa,ya) = P{x,y) 

for all x,y G E^ , then tt is exchangeable, that is, for every permutation a of 
{1, 2, . . . , N}, we have n{xa-) = 7r(a;) for all x <E E^ . 

Proof. Given a permutation a of {1, 2, . . . , A^}, define the distribution tt^ on 
E^ by TT^{x) := tt{x^). Then 

T^aiy) = ^i^)Pi^^y<y) = X! T^{Xa)P{Xa,ya) ^ ^ TVa{x)P{x,y) 

xG-E" xG-E" xG-B" 

for all y G E^ , hence by the uniqueness of stationary distributions, tt^ = tt. □ 
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Table 1: One-step transitions using the second approach, for both game A' and 
game B. From state (riQ, rii, 712), a transition is made to state (rip, n'l, n^). 



(no - 2,ni + l,n2 + 1) 
{no - l,ni - l,n2 + 2) 

(rio,ni,ri2) 

(no,ni,ri2) 
(no + l,ni - 2,712 + 1) 
(no + 2,ni - l,n2 - 1) 
(no - l,ni + 2,n2 - 1) 

(?io,7ii,n2) 
(no + l,ni + l,n2 - 2) 

(no - l,ni + l,n2) 
(no - l,ni,n2 + 1) 
(no,ni l,n2 + 1) 
(no + l,ni - l,n2) 
(no + l,ni,n2 - 1) 
(no,ni + l,n2 - 1) 



type of 

type of game winner 
player played / result 






A' 








A' 


1 





A' 


2 


1 


A' 





1 


A' 


1 


1 


A' 


2 


2 


A' 





2 


A' 


1 


2 


A' 


2 





B 


win 





B 


lose 


1 


B 


win 


1 


B 


lose 


2 


B 


win 


2 


B 


lose 



probability 



[iV(7V - l)]-ino(no - 1) 
[A^(A^- l)]-inoni 
[A^(A^- l)]-inon2 
[7V(A^- l)]-inino 

[N{N - l)]'^ni{ni - 1) 
[A^(A^-l)]-inin2 
[A^(A^-l)]^in2no 
[A^(A^-l)]-in2ni 

[A^(A^-l)]-in2(n2-l) 

N-^mqi 
N~^n2P2 
N~^n2q2 



We would like to apply results of Ethier and Lee (2009) to game B and to 
the mixed game. (They do not apply to game A' because the one-step transition 
matrix P^'^ is not irreducible, but the behavior of the system is clear in this 
case.) We restate those results here for convenience. 

Consider an irreducible aperiodic Markov chain {Xn}n>o with finite state 
space S. It evolves according to the one-step transition matrix P = {Pij)ij^^. 
Let us denote its unique stationary distribution by tt = (7ri)igs . Let w : E x E 1— )• 
R be an arbitrary function, which we write as a matrix W = {w{i^ j))ijizs and 
refer to as the payoff matrix. Finally, define the sequences {^n}n>i and {S'„}„>i 

by 

^„ := ■u;(X„„i,X„), n>l, (4) 

and 

Sn:^^i + ---+^n, n>l. (5) 

Let n denote the square matrix each of whose rows is tt, and let Z := {I — {P — 
n))^^ denote the fundamental matrix. Denote by P (resp., P) the Hadamard 
(entrywise) product P o W (resp., P o W o W), and let 1 := (1, 1, ... , 1)^. 
Then define 

/i := nPl and := tvPI - {nPlf + 2ttP{Z - H) PI. (6) 
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Theorem 2 (Ethier and Lee 2009). Under the above assumptions, and with the 
distribution of Xq arbitrary, liin„_>.oo fT-""'^E[S'„] = fi, 

> fi a.s., 

n 



limTi_j.oo ?^ ^Var(5„) — , and, if a"^ > 0, 

Sn - nfi 



dN{0,l). 



If II ~ Q and > 0, then — oo = liminf„_yoo Sn < limsup„_).Q^ S'„ = oo a.s. 

We apply this result first with S := S^r and P := Pj^\ which is clearly 
irreducible and aperiodic. We claim that the stationary distribution tt^^ is the 
iV-fold product measure tt x tt x • • • x tt, where tt = (ttq, tti, 7r2) denotes the 
stationary distribution of the three-state chain in Si with one-step transition 
matrix 

po go ' 
= h?i pi 



.P2 92 



Indeed, 

X 

N 



■ X! '^x,PB^\iyl,■■■,y^-l,X^,y^+l,...,yN),y) 



1=1 



N 



i=l 

where the first equality holds because state y can be reached in one step only 
from states x that differ from y at exactly one coordinate. Alternatively, we 
could take S := Ejv and P := .Pg^'' from Table [T] In this case the unique 
stationary distribution is multinomial(A'^, tt). 

Next, let us determine the value of /.i in the theorem. We have 



Mb 



\no,ni,n2j 

(no,ni,«2) -1-/12(^2-92) 

= 7ro(po - go) + 7ri(pi - gi) + 7r2(p2 - 52) = /-i^^ = 



N 
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because the parameterization ^ with e = was chosen to ensure the last 
equaUty. 

Now we apply the theorem with S ;= Sat and P := = "jP^^ + 

(1 — ^)P^\ where < 7 < 1, which is also irreducible and aperiodic (because 
Pj^'' is). Here the unique stationary distribution 7r|;^|_^^ is complicated. For 
example, in the simplest case, 7 = 1/2 and N ^2, 

^(1/2,1/2) 0) = (1 + + 47/9 + 60p^ + ATp^ + ilp^)/d, 

^V2.i/2)(0' 1) = ^V24/2)(l' 0) = 2(1 + p){l + p2)(ii + 15p + 9p2 + I9p3)/d, 
'^(i/2,i/2)(0' 2) = ^V2.i/2)(2' 0) = 2(1 + P)(l + P')(19 + 9/7 + 15/72 + llp^)/d, 
^(V2,i/2)(1' !) = (! + + 21/3 + 48p2 + 5^ + 27p* + 42/?5)M 
'^(V2,i/2)(1' 2) = ^V2.i/2)(2. 1) = 6(1 + P?{1 + p')(4 + P + ^p')ld, 
'^(V2,i/2)(2' 2) = (1 + p)(42 + 27p + 59/j2 + 48^^ + 2lp^ + 19p5)/d, 



where d := 2(13 - 2p + 13p2)(10 + 20p + 21/)^ + 2Qp^ + lOp"^). In particular, 

,(2) 

(1/2,1/2) 



f2) 

each entry of tt.. ,„ . is the ratio of two degree-6 polynomials in p. In another 



simple case, 7 = 1/2 and = 3, each entry of 1/2) ratio of two 

degree- 14 polynomials in p. Fortunately, explicit formulas such as these are 
unnecessary to evaluate p]^\_^y 

Let 7r|;^|_^j denote the corresponding stationary distribution on Ejv. Then 
the mean profit per turn to the ensemble of players is 

{N) ^ ^{N) p{N) 
'^(7,1-7) (7,1-7) (7,1-7) 

AT 

= (1 - 7) E ^(^!l7) (^1 , • ■ • , E ^^'^p^^ - 1=^^ ) (7) 

X i=l 

= A^"^(l-7) E ^[y\-y)in-o,ni,n2)[no{po-qo)+ni{pi-qi) 
("o,ni,«2) +^2(^2-92)] 
= A^"^(l - 7){E_(«) K](Po - 9o) +E_(]v) [ni]{pi-qi) 

(-y,l — -y) (-y,l — -y) 

-l-E (N) [?^2](P2 - 92)}- 

Now by Table [U we can compute 

, -2no{no - 1) - noni + ni{ni - I) + 2nin2 - n2nQ + n2{n2 - I) 
Efrio - no] = 7 



+ (1-7)- 



N{N - 1) 
-nopo - noqo + niqi + ^2^2 



N 

7(A^ - 3no) + (1 - 7)[no(-l) + niqi + ?i2P2] 
N 
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Similarly, 



E[n[ — rii 



_ 7(iV - 3ni) + (1 - 7)["-oPo + + '^292] 

EK - 712] = 



N 

j{N - 3712) + (1 - 7)["-ogo + niPi + "-2(-l)] 



N 

In each of these equations, we have used uq + ni + n2 = N to simplify, with 
the result that all the quadratie terms caneel and the right sides are linear in 
(no, rii, 712), at least if we replace the N in the numerators by no + ^1 + '^2- 
Next we take expectations with respect to 7r|^|_^-| to obtain 

(0,0,0) = (E_(«) [no],E_(«) [ni],E_(«) K] 

(-7,1-7) (7,1-7') (T'-I-T') 





-2 


1 




H 




-2 


:) 
















Pa 


(1- 


7) 




-1 






V P2 


92 



-1^ 

which with E_(N) [no]+E_(jv) [?7i]+E_(jv) [n2] = iV uniquely determines 

(7,1-7) (7,1-7) (7,1-7) 

the vector (E_(jv) [no],E_(iv) [ni],E_(]v) [n2]) because the matrix within 

(7,1—7) (7,1-7) (7,1-7) 

brackets is an irreducible infinitesimal matrix. Substituting into ([7]) and using 

our parametrization ([2]) with e = 0, we obtain 
(N) ^ 37(1~7)(1"P)'(1 + P) 

^^(7,1-7) 2{l + p + p2)2 + ^(5 + lOp + 6p2 + 10p3 + 5p4) + 272(1 +p + p2)2 ' 

(8) 

which does not depend on N and is positive if < p < 1, zero if p = 1, and 
negative if p > 1, indicating that the Parrondo effect is present, regardless of 
7 e (0,1), if P 7^ 1- (In the case p > 1, the effect is sometimes referred to 
as a reverse Parrondo effect. We will not make this distinction.) Temporarily 
denoting p^^^\_^-^ by p|^|_^)(p) to emphasize its dependence on p, we note that 

a fact that can also be proved probabilistically (Ethier and Lee 2009). 
When 7 = 1/2, this reduces to 

(A^) ^ 3(l-p)3(l + p) 

M(i/2,i/2) 2(10 + 20p + 21p2 + 20p3 + 10/) • 

As we will see in Section [3 this formula appears elsewhere in the literature of 
Parrondo's paradox. 



3 An alternative approach 

The method used in Section [2] to find does not extend to finding the 

variance i<j')^\_ n)2. However, a method that does extend is based on the 



8 



observation that the components of the iV-dimensional Markov chain controlhng 
the mixed game are themselves Markovian. 

For example, when game B is played, the Markov chain for player i (one of 
the N players) has one-step transition matrix 



{N-l)l3]. 



(9) 



On the other hand, the redistribution game A' affects player i only if i is chosen 
as the donor or as the beneficiary (probability {N — l)/[iV(iV — 1)] = 1/N for 
each). This leads to 



(10) 



where P^"^ denotes the one-step transition matrix for the original one-player 
Parrondo game A (not A'). In both displayed matrices, the superscript il,N) 
is intended to indicate that the underlying Markov chain controls one of the N 
players. 

From these one-step transition matrices we calculate 



A ' 



and 
With 



pil-N) ^-1 p(l) 



B ' 



p 


- IP^' - 


Ki 


-i)P'i 


N) 


p 


-ipi'r^^^ 


Ki 


~i)P'i 


N) 


p 


-IP^'^ 


Ki 




N) 



and with tt, IT, and Z chosen accordingly and 1 := (1, 1, 1)^, we have 



(7,1-7) 

The mean is readily evaluated to give 



(11) 



37(1-7)(1-P)'(1 + P) 



2(1 + p + p^Y + 7(5 + lOp + 6p2 + 10p3 + 5p4) + 272(1 +P + p2)2 ' 

which is consistent with ^ and docs not depend on N . The variance ((t|^'^^-))^ 
is also easily evaluated but is complicated; instead we provide its asymptotic 
value as iV — )■ 00 (cln ~ bjs; if limAr_j.oo o-n li>N = 1): 



9[8(1 + 7V(1 + P + P'r 
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+ 4(7 + 7^)(1 + P + + 2p + p2 + 2p3 + p'i)(l + 2p+ 12p2 + 2p^ + p^) 

+ 6(7^ + 7^)(1 + p + p2)2(3 + 20p + 30p2 + 40p3 + 66p^ + 40^^ + SOp^ 

+ 20p^ + 3p«) 

+ (7^ + 7"*) (59 + 306p + 864p2 + 1738p^ + 2781^"* + 3636p^ + 3912p^ 
+ 3636p^ + 2781p* + 1738p^ + 864pi" + 306p" + 59^^^)] 
/{A^[2(l + 72)(1 + p + p2)2 + + + 6p2 + lOp^^ + 5p4)]3}. (12) 

4 Variance parameter for game B 

Let P be the one-step transition matrix for an irreducible aperiodic Markov 
chain, let tt be its unique stationary distribution, and let 11 be the square matrix 
each of whose rows is tt. Denote by Zp := {I — {P — n))~^ the fundamental 
matrix of P. 

Lemma 3. For each positive integer N , Z(^i/isi)p+(i-i/N)i — n = N{Zp — 11). 

Proof. The one-step transition matrix {1/N)P + (1 — ^/N)! has the same sta- 
tionary distribution tt. hence the same II, so 

Z^^,N)P+ii-i,N)i = {I-[{l/N)P+{l-l/N)I-U])-^ =N{I-{P-NU))-\ 
hence it suffices to prove that 

(/ - (p - Nii))-^ - {i/N)u = {i -{P- u))-^ - n. 

For this it is enough that 

{I -{P- NU))[{I - (P - NU))-^ - (i/iv)n] 

= (/ - (P - n) + (TV - i)n)[(/ - (P - n))-i - n] 

or 

/- (i/A^)(/- (p-ivn))n 

= /-(/- (p - n))n + {N- i)n[(/ - (p - n))-^ - n]. (13) 

Now UP = pn = n, n2 = n, and so n = n(7 (p n)) and n(/ (p - 

n))"i = n. So dni) is equivalent to 

I - (i/iv)(n - (n - Nn)) = / - (n - (n - n)) + (iv - i)(n - n) 

or 7 — n = 7 — n, hence (fT5|) . and therefore the lemma, is established. □ 

We want to use this to evaluate the variance parameter for Toral's A^-player 
game B, in which there is no redistribution of wealth. The state space is Sat 
and the one-step transition probabilities are as previously specified. We assume 
the parameterization ([2|) with e = 0. 
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We have seen that the stationary distribution tt^^ is the A''-fold product 
measure tt x tt x • • • x tt, where tt = (ttq, tti, 7r2) denotes the stationary distribu- 
tion of the three-state chain with one-step transition matrix Pg"^ ■ Specifically, 

^" 2(l+p + p2)' 2(l+p-fp2)' ^2 2(l+p + p2)- 

In principle, we could use the formula :— ttPI — {ttPI)"^ + 2ttP{Z — n)Pl, 
but the evaluation of the 3^ x 3^ fundamental matrix Z is difficult, so we take 
a different approach. 

The key observation is that each coordinate of the A^-dimensional Markov 
chain is a onc-dimensional Markov chain with one-step transition matrix (|9|) or 



Pi,^-^' := {l/N)Pj^^ + (1 - 1/N)h. 



Further, the coordinate processes are independent if their initial states are, and 
they are if the initial state of the A^-dimensional process has the stationary 
distribution tt^-* on Ejv- 
As already noted in Section^ P^J'^'> = {l/N)P^J'> andP^^'^' = (l/iV)P^^'. 
By LemmaEl Z^g'^^ n = N{Z^g^ - H), so (since = iV^Vs^ = 0) 

(^^■^V := 7rP(^'^)l + 27rP(i'^)(Zl,i-^) - n)P^^'''h 

= N-^[TTP^gh + 2'KP^g\z^g^ - n)p^^^i] 

Finally, let Sn denote the profit to the ensemble of N players after n plays of 
game B, with s\i denoting the profit to player i. Then Sn = 
and the summands are independent (assuming the stationary initial distribution 
mentioned above), hence 

(cj^f) )2 ^ lini n-^Y&r{Sn) = N lim n-iVar(S'W) 



71— >-00 

2 



yielding a simple and explicit formula for ((t^'')2, which does not depend on N . 

5 Variance parameter for random mixtures 

With Sn denoting the profit to the ensemble of N players after n plays of the 
mixed game, let denote the profit to player i (one of the N players) after n 
plays of the mixed game. Then 

N 

Sn = S]j!^ , 

i=l 
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so 

Var(5„)-^Var(5W) + 2 ^ Cov(5W,5W) 

1=1 l<i<j<N 

= 7VVar(5W) + 7V(7V - l)Cov(5W , ^jfl). 
Dividing by n and letting n — )> oo, wc find that 

(.j-L.))^ = + ^(^ - (15) 

where the last superscript is intended to indicate that the underlying Markov 

chain controls players 1 and 2 of the N players. We know how to evaluate 

(it) , ' J , SO it remains to hnd cr, , ' . 
^ (7,1-7)'' ' (7,1-7) 

For this we will need an extension of (H])-®. With the same assumptions 
on {Xn}n>o (an irreducible, aperiodic, finite Markov chain in S with one-step 
transition matrix P and unique stationary distribution tt), we let w^^\'w^'^^ : 
E X E n- i? be two functions with W^^^ and W^'^^ denoting the corresponding 
matrices, and define 

:= «;W(X„_i,X„), ^[2](x„_i,X„), n>l, 

and 

5W:-Cr'+---+4', 4^1:=er'+---+d^ n>l. 
Let n and Z be associated with P in the usual way. Denote by Pl^l, Pl^l, and 
p[i'2l the Hadamard products PoW^'^\ PoW^'^\ and Po Wl^l o W^'^\ resp., 
and let 1 (1, 1, . . . , 1)"^. Then define the covariance parameter 

^[1,2] :^7rP[i'2]i_(7rp[ill)(7rP[2ll) 

+ Trpi^i (z - n)p[2i 1 + 7rp[2i (z - n)p[ii 1. 

The interpretation of this parameter is as follows. 

Theorem 4. Under the above assumptions, and with the distribution of Xq 
arbitrary, 

hm n-iCov(5W,5[2l)=a[i^2]^ 

Proof. The proof is similar to the proof that lim„_).oo 7T.^^Var(S'„) = in The- 
orem [21 which is just the special case w'^l = w'^l = w. □ 

Wc now want to apply this to find cr^^'^}^^^^ . This involves only players 1 
and 2, for which we need only a (9-state) Markov chain in E2. The reduced 
model that does not distinguish between the players but only counts how many 
players of each type there are is insufficient. 

Thinking of (i, j) € E2 as the base-3 representation of the integer 3i + j, we 
order the elements of E2 by their values (0-8). The one-step transition matrix 
for the profit to players 1 and 2 when N players are playing game B is 

P^^^^^ 7V-1[2P1'' + iN- 2)/9], 
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where ' is as in Section [2] with N = 2. The superscript (2, N) is intended to 
indicate that the underlying Markov chain controls two of the N players. The 
one-step transition matrix for the profit to players 1 and 2 when N players are 
playing game A' is 

pj2,JV) [^(^ „ 1)]-^[2Pao + 4(iV - 2)Pa, +{N-2)(N- 3)Ig], 

where Pa^ is a 9 x 9 matrix with two entries (each equal to 1/2) in each row, 
corresponding to one-unit transfers 1 — > 2 and 2^1; similarly, Pa^ is a 9 x 9 
matrix with four entries (each equal to 1/4) in each row, corresponding to one- 
unit transfers 1— >1,2— and • — >■ 2, where • represents the players 
other than 1 and 2. The functions w^^I and iijl^l can be specified as follows. 
Corresponding to matrices Pg"^ and Pai, the function w^^^ is 1 at (1 wins) and 
at • — 1; it is —1 at (1 loses) and at 1 •; and it is at (2 wins) or (2 loses) 
and at • 2 and 2 — > •. Corresponding to matrix Pao, the function lyl^l is 1 at 
2 — >■ 1; it is — 1 at 1 — > 2. The function wl^l is defined exactly in the same way 
but with the roles of 1 and 2 reversed. 

From these one-step transition matrices we calculate 

(p^z.A'))!!] ,^ 2N~\pj^^)^'\ (P^2'^))[21 := 2N-\PPf\ 

(p(2,Ar)^[l] _ [^(^_i)]-i[2(p^j[i]+4(7V_2)(P^jW], 
(P(^'^))[21 := [N{N-l)]-^[2{PAj'^+4{N-2)(PA,fl 

(p(2,A'))[i.2] :=2[7V(iV-l)]-i(P^j[i'2l. 

With 

P:^7Pi^^"^ + (l-7)Pf"\ 
pW :=^(p(^■^))[ll + (l_^)(pO.■^'))W, 

pPl :=^(p(^-^))[2l + (i_^)(p(;..A^))[2l^ 

p[1.2] ,^ ^(p(^'^))[1^2] ^ _ ^^(p(2,Ar)^[i,2]^ 

and with tt, II, and .Z chosen accordingly and 1 (1, 1, 1)^, we can evaluate 
ajlifl'^) 7rP[i^2];^ _ (ttPIH l)(7rP[2l 1) 

+ 7rpW(z - n)p[2ii + 7rp[2i(z - n)pWi 

as a function of N, at least if we fix p and 7. 
With p = 1/3 and 7 = 1/2, we conclude that 

(<^(i/2,i/2))^ = 27(-36821493886409-f 717242606475537V- 46282959184439iV2 
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+ 9902542819695iV3) (16) 
/[8331019058(-269171 + 5243477V - 3383817V2 + 72A05N% 

which is monotonically increasing in > 2, ranging from 



, (2) ,o 114315959583 „ , 

cr 1 = « 0.442636 

V (1/2,1/2)'' 258261590798 



*° r ^ (N) ^2 5941525691817 „ 

(-(i/li/2))^ = 13404609664322 " ^'^^^^^^^ 

Let us summarize our results for random mixtures. Let Sn be the cumulative 
profit after n turns to the ensemble of iV > 2 players playing the mixed game 
^A' + (1 — j)B, where < 7 < 1. Wc assume the parameterization ([2]) with 
e = 0. 

Theorem 5. If j = 1 so that game A' is always played, then P(S'„ ~ for all 
n > 1) = 1. 

= so that game B is always played, then {Sn — '5'ri-i}Ti>i satisfies the 
SLLN and the CLT with mean and variance parameters /i^' = and (cr^^)^ 



as in 

// < 7 < 1 so that both games are played, then {Sn — '5'ri-i}n>i satisfies 
the SLLN and the CLT with mean and variance parameters as in 

(or Ul\) ) and at least when p = 1/3 and ^ = 1/2, as in M0i) . When 

P7^1/3 or 77^1/2, we implicitly assume that {<y[^l_^)? > 0. 

Proof. The first conclusion is obvious. The second and third conclusions follow 
from Theorem 1, though the mean and variance parameters are obtained not 
from the theorem but by using the methods described in the text. □ 

To compare our results with those of Toral (2002), we must restore the bias 
parameter e > 0. For simplicity, let us take 7 = 1/2, as he did. Then 

^"(^2,1/2) {^[2(1 - p)^(l + p) - e(13 + 26p + 30p2 + 26p^ + 13/) (17) 
+ e'-{l - pf{l + p)- 2e-'{l + pf{l + p2)]}/{2[2(10 + 20p 

+ 21p2 + 20p^ + lOp^) - £(1 - pf{l +p)+ 3£2(1 + p)2(l + p2)]|_ 

Toral reported a simulation with p = 1/3, 7 = 1/2, e = 1/100, and N ~ 200. 
Actually, e = 1/1000 was intended (personal communication 2011). With p = 
1/3 and £ = 1/1000, ^ reduces to 193387599/6704101000 « 0.028846, with 
which Toral's estimate, 0.029, is consistent. 



6 Mean profit for nonrandom patterns 

Toral (2002) omitted discussion of the case in which his games A' and B are 
played in a nonrandom periodic pattern such as A'BBA'BBA'BB ■ ■ ■ . Let us 
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denote by [r, s] the pattern {A'YB'^ repeated ad infinitum. We would like to 
apply the results of Ethier and Lee (2009) to the pattern [r,s], showing that 
the Parrondo effect is present for all r, s > 1. (Unlike in the original onc-player 
Parrondo games, the case r = s = 1 is included.) We do this by showing that the 
mean profit per turn for the ensemble of players, Mj^jj, is positive if < p < 1, 
zero if p = 1, and negative if p > 1, for all r, s > 1 and N > 2. As we will see, 
here the mean parameter depends on N and it takes a particularly simple form 
in the limit as TV — )■ oo. 

First, Theorem 6 of Ethier and Lee (2009) is applicable. (The assumption 
there that Pa is irreducible and aperiodic is unnecessary.) But again it is 
simplest to apply the results to one or two players at a time, as we did in 
Sections [3] and [5l Let us begin by finding the mean parameter A'ljJ^ij- 

For the original one-player Parrondo games, in which 

1-1 
-1 1 
^1-10 

Ethier and Lee (2009) showed that 

1 / 1 — 1 — 

^J■[r.s] = — : — T^s.rR diag s, LC, 

r + s \ l-eil-e2/ 

where TTs^r is the unique stationary distribution of PgP^, R is the matrix of 
right eigenvectors of Pb, ei and 62 are the nonunit eigenvalues of Pb, L := 
, and C [Pb ° W)l- They further showed that this formula reduces 
algebraically to 

where 

Er,s ■■= 3a4[2 + (3a,. - l)(e"i + - 2e?e^) - (e^i + e^)](l - p)(l + p)S 

+ - eD[5(l + pfil + p') - V]}(1 - pf (18) 

and 

Dr.s ■■= 4(r + s)[l + (3a,. - l)e^][l + (3a,. - l)e^](l + p + p^fS (19) 



with ar := (1 - (-l/2)'~)/3 and S := y/ {I + p^){l + Ap + p"^). 
We apply these results but with Pa and Pb replaced by 

P^!^^) N-\2P^2^ + (iV - 2)h] and P'^'^'^ N'^P'}^'^ + [N - 

Now (P^i'^^^Y is given by the same formula as P\ but with redefined as 

a,. [1 - (1 - 3/iV)n/3, (20) 
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and (P^^'^-*)* has the same spectral representation as Pg but with the nonunit 
eigenvalues replaced by 

er:^l-i^, e. :^ 1 - i^, (21) 
where e J and 62 are the nonunit eigenvalues of Pb , namely 



' ■ 2 2(l + p)(l + p2)' ^ • 2 2(l + p)(l + p2)- 

The matrices R and L arc unchanged. 
We conclude that 

= NEr.s/Dr,s, (22) 

where E^^s and Dj.^s are as in and (|19p with only the changes (PU)) and (PT|) . 
For example, this leads to 

= 3N{2N - 3)(1 - pf{l + p)/{2[18(l + p + p^f - 3N{13 + 26p 

+ 30p2 + 26p3 + 13p*) + 2N^{10 + 20p + 2lp^ + 20^^ + lOp"*)]} 

and 

^[1,2] 

= 27V(1 - pf{l + p)[-3(l + p + p^f + iV(10 + 20p + 21p2 + 20p3 + IQp'^) 

- 97V2(1 + p)2(l + p2) + 3^3(^ ^ + p2)]/[36(l + p + /92)4 

- 12iV(l + p + p2)^(ll + 22p + 2V + 22p3 + lip'') + A^^(193 + 772p 
+ 1660p2 + 2548p3 + 2938^"* + 2548p^ + 1660p^ + 772p^ + 193p^) 

- 3iV^(l + pf{Ai + 86p + 145p^ + 172p3 + 145p'' + 86p'^ + 43p'^) 
+ N-^{1 + p)^(35 + 70p + 113p2 + 140p'^ + 113p^ + 70p''^ + 35p^)]. 

Both of these functions arc positive for all > 2, as can be seen by expanding 
numerators and denominators in powers of A^— 2 and noticing that all coefficients 
are polynomials in p with only positive coefficients. 

Although these formulas become increasingly complicated as r and s in- 
crease, their limits as A'^ — >■ 00 have a very simple form. To see this, it suffices 
to note that 

and similarly for e| , so p2|) converges as — > 00 to 

3rs(l-p)3(l + p) 

9r2(l + p)2(l + p2) _^ 9^s(l + ^)2(l + ^2) + 2s2(l + p + p2)2 ' 

which coincides with (jH)) (or ([lip ) when 7 = r/(r + s). This limit is positive if 
< p < 1, zero if p = 1, and negative if p > 1, so we conclude that the Parrondo 
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effect is present for all r, s > 1 , as long as N is large enough and p ^ 1- This 
relationship between the random-mixture case and the nonrandom-pattern case 
is not present in the original one-player Parrondo games except in a single case 
(r = 2, s = 1). (We have confirmed this for r,s > 1 and ?' + s < 75 and expect 
that it is true generally.) 

We now verify that the Parrondo effect is always present. We begin with a 
lemma. 

Lemma 6. If < a < b < c, then (c" — b"')/{b" — a") is increasing in n > 1. 

Proof. Divide both numerator and denominator by 6" to see that wc can, with- 
out loss of generality, assume that 6=1. So the aim is to show that 

c" - 1 c"+i - 1 

< -; TT: n> 1, 

1 - a" 1 - a"+i ~ 

or that 

c" - 1 a" - 1 
— n 7 < — n 7i n > 1. 

For this it is enough to fix n > 1 and show that the function 

- 1 



- 1 ' 



defined by continuity at a; = 1, is decreasing on (0, oo). Its derivative has the 
same sign as 

-[a;"+i-(n + l).T + n], 

so it is enough that the quantity within brackets is positive for x > 1 and 
< a; < 1. First suppose that a; > 1. Then 

a;"+i - (71 + l)x + n = {x-l + 1)"+^ - {n + 1)(.t - 1) - 1 

= (-ir- + (";')(-ir + -- + (:!;)(-i)^ 

> 0. 

Next suppose that < a; < 1. Then 

x"+^ - (n + l)a: + n = x''+^ - 1 - (n + l)(a: - 1) 

= (x - l)(a-" + x"-^ + --- + X + 1)- {n + l){x - 1) 
= {x- l)[x" + + ■ ■ • + a; + 1 - (n + 1)] 
> 0. 

This completes the proof. □ 

Theorem 7. /i[^J] is positive i/0 < p < 1, zero if p = 1, and negative if p > 1, 
for all r, s > 1 and N > 2. 
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Proof. Denoting temporarily by p-^^^j^ (p) to emphasize its dependence on 
p, it can be shown algebraically or probabilistically that 

^S(i/p) = -^S(p), 

so it will suffice to treat the case < p < 1. First, |3a,- — 1| < 1 and ei, 62 G (0, 1), 
so Dr.s > 0. Since > 0, it suffices to show that 

[2 + (3a. - l){el + - 2ele^) - (e? + e^)](l - p)(l + p)S 
+ arie'2 - el)[5{l + pf{l + p^) - V] > 0. 

Discarding the — term (since e| — < 0), it is enough to show that 

(1 - e-J)[l + (3a. - l)e5] + (1 - e^)[l + (3a. - l)e\] 

.„,,_.,!<l±£)a-!)>o. ,3, 

NoweJ = (-l+x)/2 ande^ = (-l-x)/2, where a; {I- p)S + p){l + p^)] e 
(0, 1), so ei = (27V - 3 + x)/{2N) and 63 = (2iV - 3 - x)/{2N). 

Let us ffist assume that iV > 3. Then 3a,. — 1 < 0, so, replacing e\ and e| 
within brackets in ()23p by 1, we need only show that 



or that 



3(l-eD + 3(l-e^)>(e?-e^)^(i±^fc^, 

2(2iV)- - [(2jV - 3 + x)- + (2jV - 3 ~ xf] 5 

[(2A^-3 + a;)''-(27V-3-a;)'']/x 3' ^ ' 

The denominator is a polynomial of degree s — 1 in x with positive coefficients, 
while the term within brackets in the numerator is a polynomial of degree s in 
X with positive coefficients. So the left side of (|24|) is decreasing in x, and it 
suffices to verify it at x = 1. For this we notice that 

2{2NY - [{2N - 2Y + {2N - AY] _ ^ N"" - {N - lY ^ 
(2N - 2Y - {2N ~ AY " (N -\Y-{T^ -2Y ^ ' 

and the fraction on the right is increasing in s > 1 by Lemma [S] At s = 1 the 
value is 3, so the desired inequality holds. 

It remains only to consider the case N = 2. The same argument works if 
r is even because then 3a. — 1 < still holds. If r is odd, we can replace the 
quantities within brackets in ((23)) by 1 and can replace a. in the second line of 
by ai = 1/2. Thus, we need only verify (|24l) with 5/3 replaced by 5/2, and 
of course it still holds. □ 



7 Remark on a "coincidence" 

Wc can prove algebraically that 

lim u^^'^ ~ u^^^ -{■M2\a^^^ - CM 2^a^^^ - u^^' 

M-5oo'^['"^''l ~ ^(1/2,1/2) ~ '>"'/^/^^(2/3,l/3) ~ ^'^ I ^ ' [2 ,1] ~ 
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^ 3(l-p)^(l + p) 

2(10 + 20p + 21p2 + 20p3 + 10p4) 

for all r > 1 and N > 2, where superscripts refer to the number of players. 
(For superscripts equal to 1, the games are A and B, the original one-player 
Parrondo games. For superscripts 2 or larger, the games are A' and B.) The 
first equality is from Section [6l Can the others be explained probabilistically? 
We can elucidate at least the second equality, while the third and fourth remain 
partially unexplained. 

Since does not depend on N, it is enough to verify the identity with 

N = 2. Let us consider the profit of one player when two players are playing. 
Recalling Q and with iV = 2, we have 

PiP):=Pl^) and P^^) := (l/2)(p(^) + /a). 

The former is just the one-step transition matrix for the original one-player 
game A, and we have 

(l/2)PiV^) + (1/2)P(^'^) = (1/2)p1'^ + (1/4)p4'^ + (1/4)J3. 
The left side describes the (1/2,1/2) mixture of games A' and B, as viewed 

(2,) 

by one of two players. Its mean is {1/2) fi^^^'^^ ^^^y The right side describes the 
(2/3, l/3)-mixture of games A and B if we ignore the (l/4)/3 term and normalize 
to ensure a stochastic matrix. That term just slows down the process, making 
one- fourth of its transitions null. So its mean is (3/4)/i|2/3 1/3)- These arc equal, 

so A^|]^^2 1/2) ~ (^/2)/^(2/3 1/3)' claimed. 

This can be regarded as a more correct version of the argument sketched in 
the third paragraph of page L307 of Toral (2002) and attributed to an anony- 
mous referee of that paper. 



8 Variance parameter for nonrandom patterns 

Wc can evaluate the variance parameter (cj^^j)^ for the nonrandom pattern 
[r, s] in the iV-player games directly for small N , using the state space T,n with 
its (^+^) states. We apply (25)-(27) of Ethier and Lee (2009), obtaining, for 
example, 

[1,1]'' 

= [9(466 + 2680p + 7621^^ + 16310/)^ + 29018^" 4- 41582p^ + 51471p^ 
+ 55998p^ + 51471p*^ + 41582^^ + 29018^^° + 16310p" + 7621^^2 
+ 2680p" + 466p^^)]/[4(2 -p + 2p^){l0 + 20p + 21p^ + 20p^ + Wp^f], 

which reduces when p = 1/3 to 74176355601/141627323986 « 0.523743. Since 
N ~ 2, this is a computation involving 6x6 matrices. 
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To get results for larger N, we apply the method of considering one or two 
players at a time. By analogy with (jl5|) . we have 



where 



u=0 

s-1 



t;=0 

2 J2 ^pxPAipr^-' -npx)PAi 



0<.u<Cv<r — l 
r—1 s—1 



+ 2 5] 7rPlPi^Ps(Pr"-'-nPiP^)Psl 



0<ii<i;<s-l 
- r—1 7^ — 1 

2 

r—1 s—1 



^ ^ 7rPiP^p;-"-ip|(z - n)pj(p^i 

r—1 s—1 

J2 E ^^'A^'APr-'^PK^ - n)p;pi^PBi 

u=0 t>=0 
s-1 r-1 

^ ^ nP^^Pi^PsPr^-'iz - n)PXPAi 

11=0 D=0 

s-1 s-1 

^ ^ 7rP;Pi^PBPr""'(^ - n)PlPi^PBl 



(from Ethier and Lee 2009) with Pa and Pg replaced by P^/^-* and P^^'^-* as 
defined in Section [3l 

The covariance term, requires an extension of the preceding for- 

mula to covariances. We omit the details of the derivation and just give the 
result: 

([1,2], iV) _ ■ ' ' ^ 



(T, 

r,s 



-^\j2[nPXP!^''h~ i7rP:iP]^h){r:PXP!^h)] 

^ * u=0 

s-1 

+ ^[TTP^Pi^Pi^'^ll - (7rP;iPi^pWl)(xP;Pi^P|Il)] 

^ [ttp^pJ^i (p;(— 1 - npi)pi2i i 

0<u<u<r-l 



■u=0 
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+ 7rPipjfi(pr""' - nP^P^i] 

r—1 s—1 

.^^[7rPipW(pr"-' -nPi)Pi;p|ii 

u=0 v=0 

+ 7TP;(pjfi(p;-"-i - np;i)Pi^pWi] 
^ [7rPiPi^pW(p;^""-i-nPiPi;)p|ii 

- - + TrPiPi^pfi (p;^-"-i - np;p^ )pw 1] 

r — 1 r — 1 

. 5] 5][7rPipWp;-"-ip|(z - n)Pipjfii 

+.PIpifipr"-PB(z-n)p;(pWi] 

X: E[^J'A^^'^'r""'J'B(^ - n)PiPi^p|ii 
"=°^=° + 7rPj(pjfip;-"-ip|(z - n)p;p|pWi] 

s-1 r-1 

■ 5]5][7rPIP^pWpr "-^(Z - U)PXP^^h 

+^p;p;^p|ip-- ^(z-n)PipWi] 

s-1 s-1 

EE[7rP;Pi^pWp-«-i(z - n)p;Pi^p|ii 

+ 7rp;pi^p|ipr"-'(z - n)p;pi;pWi] 



with Pa and Pg replaced by P^^'^-* and Pg'^^ as defined in Section [5l 
By analogy to p6|) . with p = 1/3, we conclude that 

(cr[fj)j)2 ^ 9(615639408424560- 64089266202140407V + 295415459578941397V2 

- 802148142000374917V^ + 143582273075781927A'"' 

- 1791925578025431307V5 + 160434481099881996^^** 

- 104152159483211664iV^ + 487990916854784687V*^ 

- 161375219565952467V^ + 3584898779152593iVi° 

- 4816333990183977V" + 29679648590925iVi2) 

/[2(1521 - 3174iV + 1609iV^)^ (3285360 - 98161207V + 12525387iV^ 

- 87255897V^ + 35019287^ - 7688517V^ + 724057V'*)], (25) 

which is monotonically decreasing in A^. ranging from 

, (2) 74176355601 
(--[Ml) - 141627323986 

to 

, (AT),. 5935929718185 

lim cr , V,) = « 0.442827. 

AT^oo^ [1.1]^ 13404609664322 
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This last number differs slightly from the corresponding limit in the random- 
mixture case, showing that the variances lack the nice property that the means 
enjoy. 

Similar formulas can be found for other [r,s], assuming p ~ 1/3. With 
[r,s] = [1,2] (resp., [2,1], [2,2], [2,4], [3,3], [4,2]), we get in place of ^ a ratio 
of degree-21 (resp., 24, 33, 51, 54, 57) polynomials in TV, and 



1891312136577 
6060711605323 



0.312061 if [r,s] = [2,1] or [4,2], 



N 



= mS43l - 0-442827 if [r,s] = [1,1] or [2,2] or [3,3], 



136286243910 
252688187761 



0.539346 if [r,s] = [1,2] or [2,4] 



In particular, it appears that the result for [r, s] depends on r and s only through 
r/(r + s). We have confirmed this only in the several cases shown above; a proof 
for all integers r,s > 1 seems difficult. 

Let us summarize our results for nonrandom patterns. Let Sn be the cu- 
mulative profit after n turns to the ensemble of A'^ > 2 players playing the 
nonrandom pattern (A^YB" (denoted by [r,s]) repeatedly, with r and s being 
positive integers. We assume the parameterization ([2]) with e = 0. 

Theorem 8. {Sn — 5'Tt-i}n>i satisfies the SLLN and the CLT with mean and 
variance parameters fJ-\^\ computable for specified r, s > 1 as a function of 

N > 2 and p > and (fjj^j,')^ computable for specified r, s > 1, N > 2, and 

p > 0. We implicitly assume that (cjj^^j)^ > 0. 

Proof. The proof is as for Theorem [5l except that, instead of citing Theorem 1, 
we cite Theorem 6 of Ethier and Lee (2009). □ 



9 Sample variance of players' capitals 

Recall our notation in 
after n trials, so that 



Recall our notation in which S"!*' is the capital of player i (one of the N players) 



N 

On ■— On 

is the total capital of the N players after n trials. Toral (2002) simulated the 
sample variance of , • • ■ , S^n^ , namely 

r JV , N , 2n 

L 1=1 ^ i=i 

finding that it grows approximately linearly in n. Let us replace this sample 
variance by its unbiased (at least in the case of a random sample) version. 



|- JV , N ,2 



TV- 1 



1=1 ^ 1=1 
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and let us consider its expectation E[(a-'^''^^)^] instead of the random variable 
itself. We can evaluate this using the results already obtained. Indeed, 



1 



N- 1 
1 

N- 1 
1 

N - 1 
1 

n- 



Ee[(5W)2]-1e[(5„ 

1=1 

^ 1 

E{Var(5H) + (£[5^])^ - -{Var(5„) + (E[5„])n 

1=1 

7VVar(5W)-lvar(5„) 



iV- 1 

as n — 7> oo, hence 



,([1.2], W)l 



hm n-iE[(aW)2] = (a(i'^))2 - a^Ii-^l.^). 

n— >-cxD 

We have omitted subscripts A' , B, (7,1 — 7), and [r, s] because we want to apply 
this formula in all cases. 

Let us first consider the random-mixture case with p arbitrary. With 7 = 1/2 
we have 



V''^(l/2.1 



(1/2,1/2)^ 



'(1/2,1/2) 



- 27(97 + 606/9 + 1926p2 + 4262^^ + 728^ + 989^ + 10911p^ 
+ 989V + 7284/ + 4262p^ + 1926p^" + 606/?" + 97^^^) 
/[2Ar(10 + 20/9 + 21/9^ + 20/9^ + lOp'')^] (26) 

as — 00, which is ([T2|) with 7 = 1/2. With 7 = (only game B is played) 
we have 



/„(1,JV)^2 _ „([1,2],JV) 



3p 



1 + p + p' 

Finally, with 7 = 1 (only game A' is played) we have 



■ 1 
N' 



iV- 1 



2 

N' 



As Toral found, the result for the mixed game lies between those for games 
A' and B, and this is true regardless of p > 0. (Our results arc not directly 
comparable to his because we have taken the bias parameter e to be 0.) 

Finally, let us consider the nonrandom pattern [1,1]. We find that (c^ 'jT'')^ — 



a 



([1,2],JV) 
[1,1] 



has the same asymptotic value as in (|26p. so we have another coinci- 
dence. It appears that these expected sample variances have essentially the 
same property that the means enjoy, namely that their asymptotic value for the 
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case of the nonrandom pattern [r, s] depends only on r/(r + s) and p and is equal 
to the asymptotic value for the case of the random mixture with 7 = r/(r + s). 
We have confirmed this in the six cases s > 1 and r + s < 4, but a proof for 
all integers r, s > 1 seems difficult. 
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